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ABSTRACT 


The 7 ring is one of the narrow rings of Uranus, consisting of a dense core that is 1-2 km wide and a diffuse outer 
sheet spanning about 40 km. Its dense core lies just exterior to the 3:2 Inner Lindblad Resonance of the small moon 
Cressida. We fit the 7 ring radius residuals and longitudes from a complete set of both ground-based and Voyager 
stellar and radio occultations of the Uranian rings spanning 1977-2002. We find variations in the radial position of the 
m ring that are likely generated by this resonance, and take the form of a 3-lobed structure rotating at an angular rate 
equal to the mean motion of the moon Cressida. The amplitude of these radial oscillations is 0.667+0.113 km, which is 
consistent with the expected shape due to the perturbations from Cressida. The magnitude of these variations provides 
the first measurement of the mass and density of the moon Cressida (m = 2.5 + 0.4 x 1017 kg and p = 0.86 + 0.16 
g/cm?) or, indeed, any of Uranus’ small inner moons. A better grasp of inner Uranian satellite masses will provide 
another clue to the composition, dynamical stability, and history of Uranus’ tightly packed system of small moons. 
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1. INTRODUCTION 


In March of 1977, Elliot et al. (1977), Millis et al. 
(1977), and Bhattacharyya & Kuppuswamy (1977) dis- 
covered nine narrow rings around the planet Uranus by 
measuring the light blocked by each ring before and af- 
ter Uranus occulted the star SAO158687. Since then, 
the Uranian rings have been studied extensively with 
ground based stellar occultations (Millis & Wasserman 
1978; Nicholson et al. 1978; Elliot et al. 1981b,a; Nichol- 
son et al. 1981; French et al. 1982; Sicardy et al. 1982; 
Elliot et al. 1983, 1984; French et al. 1986b,a; Elliot et al. 
1987; French et al. 1988, 1996). Occultations provide 
very precise radial locations of the rings at different lon- 
gitudes in their orbits around Uranus. French et al. 
(1988) found that the main rings of Uranus consist of 
six measurably eccentric rings (6, 5, 4, a, 6, and €) and 
three nearly circular rings (m, y, and ó). In the past, 
measurements of the 7 ring’s radius have not shown the 
ring to be anything but circular. The 7 ring also features 
a broad low optical depth sheet extending approximately 
40 km exterior to its narrow core (Elliot et al. 1983). 

During the Voyager 2 flyby of Uranus Smith et al. 
(1986) discovered ten new small inner moons, but no one 
has ever measured their masses or densities. Nine of the 
moons orbit within a radial range of 20,000 km, mak- 
ing the group one of the most tightly packed systems 
of interacting satellites in our solar system. Lissauer 
(1995) estimated the masses of the inner moons assum- 
ing densities equal to that of the larger moon Miranda 
(Jacobson et al. 1992) and shapes estimated with pho- 
tometry (Thomas et al. 1989), but stated that at least 
some of Uranus’ small inner moons are significantly less 
massive than these estimates. The lifetime of this sys- 
tem is highly sensitive to the masses of the individual 
satellites (French et al. 2015). In fact, prior to the knowl- 
edge of the even less stable moon Cupid (Showalter & 
Lissauer 2003; French & Showalter 2012), Duncan & Lis- 
sauer (1997) showed that Desdemona could collide with 
either Cressida or Juliet within the next 4 — 100 million 
years, depending on the masses of the satellites involved. 
The discovery of the dusty v and yu rings (Showalter & 
Lissauer 2006), near the orbits of Portia/Rosalind and 
Mab respectively, hints at the possibility of an evolving 
inner ring-moon system dominated by accretion (Tis- 
careno et al. 2013). Kumar et al. (2015) also argue that 
anomalies in Mab’s orbital motion may be explained by 
a ring-moon system that is undergoing re-accretion after 
a recent catastrophic disruption. 

Here we investigate a complete set of Uranian 7 
ring occultation observations spanning their discovery 
in 1977 to 2002. We find that the 7 ring’s radii ex- 
hibit a 3-lobed structure rotating around Uranus at the 


mean motion of the moon Cressida. We argue that this 
structure is a result of the 7 ring’s close proximity to 
Cressida’s 3:2 inner Lindblad resonance (ILR). One of 
the maxima in the ring’s radius aligns with Cressida, as 
expected for the stable ring structure located exterior to 
the resonant radius. The measured radial amplitude of 
this ring structure and its distance from the resonance 
allow us to estimate Cressida’s mass, and thus obtain 
the first gravity-based mass measurement of any inner 
Uranian moon. 

We have only been able to find three previous men- 
tions of the Cressida 3:2 ILR and its association with the 
n ring. Porco & Goldreich (1987) identified the most 
relevant resonances in the Uranian ring-moon system 
and made a case for Cordelia and Ophelia shepherding 
the outermost e ring through torques generated by the 
Lindblad resonances located appropriately on the ring’s 
inner and outer edges (Goldreich & Porco 1987). They 
also note single resonances that could be perturbing the 
y and ó rings. Finally, they state: “The only isolated 
first-order satellite resonances which fall near any of the 
remaining rings are located interior to the 7 ring.” Porco 
& Goldreich (1987) list both the Cressida 3:2 and the 
Cordelia 13:12 resonances, located at a = 47171.6 + 0.3 
km and a = 47173.0 + 0.3 km respectively. These reso- 
nances fall 3—5 km interior to the 7 ring. They calculate 
the widths of both resonances to be ~ 1 km and dismiss 
the possibility that either resonances is perturbing the 
n ring. Murray & Thompson (1988) later marked the 
location of the Cressida 3:2 ILR in their figure display- 
ing a radial scan of a high phase image of the Uranian 
rings acquired by Voyager 2. Subsequently, Murray & 
Thompson (1990) noted that this resonance needs to be 
re-examined using updated satellite parameters. At the 
time, with a smaller data set, there was no detection of 
either an m = 3 or an m = 13 mode in the 7 ring, nor 
any other modes due to resonances with known satel- 
lites having observed effects on any of the other pre- 
viously noted rings' edges (French et al. 1988). Thus, 
it was only sensible to dismiss these resonances, and it 
is reasonable that they have not been of interest since. 
We are only able to make this discovery now because 
we have a larger set of occultation data extending from 
1977 through 2002. 

We present the data used in this analysis in Section 
2, and describe our ring particle streamline model and 
our mode detection methods in Section 3. In Section 4, 
we report the parameters of our fit to the 7 ring and 
calculate the mass and density of Cressida. Finally in 
Section 5, we discuss potential implications for the dy- 
namical stability of the tightly packed system of inner 


Uranian moons and the possible composition of Cres- 
sida. 


2. OBSERVATIONAL DATA 


'The observational data used for this analysis consist 
of 49 individual occultation observations of the 7 ring. 
In the appendix, Table 3 contains each occultation's ring 
intercept time, inertial longitude, and mid-radius deter- 
mined using a simple square-well model for profile fit- 
ting, developed by Elliot et al. (1984) and used in later 
orbit determinations of the Uranian rings (French et al. 
1986b,a, 1988, 1991). Of these 49 observations, 46 are 
Earth-based stellar occultations, two are Voyager 2 Ra- 
dio Science Subsystem (RSS) radio occultations, and one 
is a Voyager 2 Photopolarimeter Subsystem (PPS) stel- 
lar occultation. Several of the observations are ingress 
and egress pairs from the same occultation of Uranus 
and its rings. 

For each Earth-based occultation, an instrument 
recorded the brightness of the background star as a 
function of time. As the Earth moves relative to Uranus 
the rings can block the star's light, leaving each ring’s 
mark as a sharp decrease in the recorded brightness of 
the star for some amount of time related to the width 
of the ring. Typically the observations were detected 
with an InSb photometer in the 2.2 um band, using the 
K filter, where Uranus is fainter than the rings. Most 
observations provided limited information about the ra- 
dial structure within the rings, and here we are making 
use only of the estimate of the radius of the mid-point of 
each ring occultation profile. Interested readers should 
see Elliot (1979) for a review of stellar occultation stud- 
les of the solar system and Elliot & Nicholson (1984) for 
a review of this observation method specific to the rings 
of Uranus. 

'To identify possible Uranus occultation opportunities 
Taylor (1973) compared positions of Uranus to stellar 
positions in the Smithsonian Astrophysical Observatory 
(SAO) catalog. Once the rings were discovered, it be- 
came more appropriate to utilize dimmer stars that are 
bright in the 2.2 wm band. Thus, Klemola & Marsden 
(1977) searched for stars on photographic plates con- 
taining star fields ahead of Uranus and created a list of 
ideal future occultation observations. Additional lists of 
this type were compiled by Klemola et al. (1981), Mink 
& Klemola (1982, 1985), Nicholson et al. (1988), and 
Klemola & Mink (1991). 

'The Voyager 2 PPS stellar occultation only detected 
the 7 ring on egress (Lane et al. 1986; Colwell et al. 
1990). In the case of the Voyager 2 RSS occultations, 
the RSS instrument illuminated the rings at 3.6 cm and 
13 cm wavelengths in the direction of Earth once be- 
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yond the ring plane. Stations on Earth detected the 
diffracted signal and relative phase change, to later be 
reconstructed into high-resolution radial optical depth 
profiles after the removal of diffraction effects (Tyler 
et al. 1986; Gresh et al. 1989). Presently, ground based 
occultation opportunities are rare because Uranus has 
passed out of the dense Milky Way background, dras- 
tically reducing the density of appropriate background 
stars. The rings are also no longer as open to our view 
from Earth as they were in the 1980s because the ap- 
parent aspect of the ring plane as viewed from Earth 
changes over time. 


3. RING PARTICLE STREAMLINE MODEL AND 
FITTING METHOD 


The procedure used here follows that of French et al. 
(1986b, 1988, 1991) for the Uranian rings, more recently 
employed by Hedman et al. (2010), Nicholson et al. 
(2014a,b), and French et al. (2016b) for analyses of Sat- 
urn's non-circular narrow rings, gaps, and edges. After 
taking account any inclination relative to the equatorial 
plane, the majority of narrow rings are well-fit by simple 
precessing Keplarian ellipses whose radii are described 
by: 

a(l — e?) 


"ges 1+ecos f’ 


(1) 
where the true anomaly f = À — wo — w(t — to). Here, 
the radius of the ring will vary with longitude and 
time £, where a and e are the ring's semi-major axis and 
eccentricty, wo is the ring’s longitude of periapsis at the 
time to, and c is the ring's apsidal precession rate. We 
can approximate a nearly circular (e ~ 0) ring's radii as 
r c a(1 — ecos f). 

Additionally, several rings are found to contain forced 
radial oscillations and in a few cases there are even rings 
whose structure is dominated by free normal mode oscil- 
lations. In these cases, the structures are distinct from 
circles or ellipses and their radii are described by: 


r(A,t) — a — Am cos(m6), (2) 


where 0 = À — Qp(t — to) — dm, following the formalism 
of Nicholson et al. (2014a,b) and French et al. (2016a). 
Here, the systematic radial oscillations of the rings form 
a m-—lobed figure rotating around their planet at a pat- 
tern speed Q, with a radial amplitude Am and phase 
ôm. We show some exaggerated models of m-lobed ring 
streamlines, resulting from both free normal modes and 
Lindblad resonances, in Figure 1. While individual par- 
ticles follow normal elliptical orbits, described by Equa- 
tion 1, the ring as a whole consists of streamlines with 
m azimuthally symmetric radial minima and maxima 


Figure 1. The above shapes are an exaggerated representa- 
tion of the m-lobed ring streamlines we detect in the Uranian 
rings. For each case of m, we have shown 3 streamlines with 
slightly different semi-major axes and a positive eccentricity 
gradient. Our addition of an eccentricity gradient results in 
a narrower ring width at periapsis, as is the case for several 
of the Uranian rings. 


rotating around the planet with the frequency 


(m — 1)n + Wsee 
p — 0 — a 


Q (3) 


Here, the mean motion n and apsidal precession rate 
Wesec are functions of the semi-major axis a of the ring, 
and m can be any positive or negative integer. If we 
consider the case of m = 1 we find that Qp = sec, 
A, = ae, and 6; = co, so that r is equivalent to the 
approximation of Equation 1 above. 

In the case of a free normal mode oscillation, the pat- 
tern speed will be equal to the expected pattern speed 
obtained by evaluating Equation 3 at the semi-major 
axis of the ring. However, if the ring is perturbed 
by a satellite through a first-order Lindblad resonance, 
then the ring structure will have a forced pattern speed 
matching the mean motion of the perturbing satellite 
Ns = Qp and will differ from the expected pattern speed 
slightly based on the ring’s separation from the exact 
radius of the resonance |a —a,es|. The ring is perturbed 
by the satellite due to the near commensurate ratio of 
the ring particles’ orbital periods and the period of the 
perturbing satellite. As such, first-order Lindblad res- 
onances are defined by |m| : |m — 1|, where for every 
|m| orbits of the ring particle, there are |m — 1| orbits of 
the corresponding satellite. In the majority of cases, the 
perturbing satellite lies at a larger semi-major axis than 
the ring (a, > a). The relevant resonances in this case 
are called inner Lindblad resonances (ILR) and are as- 
signed positive values of m. In the rare case of a satellite 
located interior to the rings it is possible to have both 
ILR and outer Lindblad resonances (OLR) at locations 
within the rings, allowing for negative values of m. 

'The condition for a first-order Lindblad resonance is 
that the resonant argument: 


p = m(A— As) - (A — w) (4) 


m 


is constant in time. Here À and A, refer to the lon- 
gitudes of a ring particle and the satellite respectively 


Figure 2. An exaggerated cartoon model of ring particle 
streamlines around a planet and near a 3:2 ILR with an ex- 
terior moon, in the co-rotating frame of the moon. We've 
marked the resonant radius with a dashed line and included 
three ring particle streamlines on either side of the resonance. 
This shows the stable configuration on either side of the res- 
onance, where a periapse (apoapse) is aligned with the moon 
interior (exterior) to the resonance. 


and z is the longitude of periapsis of the ring particle. 
If we consider a conjunction of the ring particle and the 
satellite (À — A, = 0) occurring when the ring particle 
is also located at its longitude of periapsis (A — w = 0), 
then the condition that y is constant implies that all 
future conjunctions will occur when the ring particle is 
near periapsis. In general, this means that the ring par- 
ticle will always be in the same phase of its orbit when 
it passes longitudinally close to the satellite. This al- 
lows the perturbing satellite to force the eccentricity and 
periapsis locations of streamlines located near the res- 
onance. In Figure 2 we show a cartoon model of the 
resulting streamlines surrounding a 3:2 ILR in the co- 
rotating frame of the perturbing satellite. Interior (ex- 
terior) to the resonant radius, marked with the dashed 
line, the streamlines are stable when oriented such that 
one of the three periapses (apoapses) is aligned with the 
satellite. 

In short, our procedure is a search for patterns in the 
varying mid-radii measurements of the rings. Each ring 
occultation observation provides the ring's radius at a 
particular longitude and time. To search for patterns in 
each ring we need the observed parameters, an m value 
to test, and the resulting expected pattern speeds for 
that m value. For each test of m, we compute the ex- 
pected patten speed for the semi-major axis of the ring 
using Equation 3 and create an array of 100, 000 pattern 
speeds, evenly spaced in increments of 0.00001? /day, 


surrounding the expected pattern speed. Using each 
pattern speed we calculate m6, for every ring observa- 
tions’ longitude À and time t, using an initial epoch time 
to of UTC 1977 MAR 10 20:00:00.00. We can then com- 
pute the observed ring radii r vs. m0 mod 360? and fit 
the data to a single sinusoid. The resulting fit parame- 
ters are a, Am, and ôm, allowing us to compute model 
values of r using Equation 2. We compute the RMS 
deviation of the observed radii and the model radii for 
each m’s 100, 000 test pattern speeds and look for a RMS 
minimum to identify the best fitting pattern speed. 

We first checked our algorithms by searching for 
known structures in the Uranian rings. In several rings 
(6,5, Jo, B, and €) we can easily detect RMS deviations 
that drop to nearly zero (sub-km) with the proper pat- 
tern speed and m input. These are the rings that largely 
follow classical Keplerian ellipses (m — 1) and whose 
pattern speeds equal the rings’ apsidal precession rate, 
Qp = Wesec. The 7,7, and 6 rings are nearly circular 
and their residuals are relatively larger when fit with 
a low amplitude m = 1 ellipse. We are also able to 
identify the known m = 2 structure of the ó ring and 
the combination of m = 0 and m = 1 for the y ring 
(French et al. 1986b). 

We decided to identify the strongest resonances in the 
Uranian rings to have a better idea of the resonantly 
forced modes that are the most likely to be detected. To 
quantify the ‘strength’ of the resonances in the system 
we chose to compare the expected forced radial ampli- 
tude on rings near each of the possible resonances in the 
main ring system. We use Equation 10.22 from Chapter 
10 of Murray & Dermott (1999), 


__2a02%(ms/my) ful 
Am = 8G = Da area ý 


where Am is the forced radial amplitude of a ring particle 
in the vicinity of a Lindblad resonance (Porco & Nichol- 
son 1987; Goldreich & Tremaine 1982). This amplitude 
is a function of the ratio of the perturbing satellite and 
central planet masses Ms/Mp, the radial separation of 
the ring and the resonance |a — ares|, the ratio of the 
ring and satellite semi-major axes o = a/as, and the 
Laplace factor fa, that depends on 7, the integer co- 
efficient of the satellites longitude in the resonant argu- 
ment, which is equivalent to m in the case of a first-order 
Lindblad resonance. As shown in Figure 10.10 of Mur- 
ray & Dermott (1999), zeli varies between 1.5 and 1.6, 
depending on j. Note that Equation 5 isn't necessarily 
applicable for all cases. If |a — ares| is smaller than the 
resonance half-width, then A,, calculated using Equa- 
tion 5 is not a good estimation of the radial amplitude 
produced by the resonance because in this regime neigh- 
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boring streamlines will cross and collisional dissipation 
cannot be ignored. 

In Figure 3 we display the forced amplitude on all 10 
rings (inner and outer edges for the < ring) due to all 
possible first-order Lindblad resonances of all Uranian 
moons out to Perdita. For the estimated mass of each 
moon, we use the radius measurements of Karkoschka 
(2001a) and Showalter & Lissauer (2006) and consider a 
range of densities from 0.5 to 1.3 g/cm?. In the left half 
of Figure 3 all resonances mentioned by Porco & Goldre- 
ich (1987) are apparent in addition to a previously unex- 
plored 2:1 ILR with Portia in the proximity of the 6 ring. 
In the right side of Figure 3 we compare the amplitudes 
of the strongest resonances over a range of moon densi- 
ties. The fainter patches in the left side of Figure 3 are 
due to resonances inducing much weaker amplitudes due 
to their large distance from the rings. Despite the sepa- 
ration in semi-major axis of the 7 ring from the Cressida 
3:2 ILR the 7 ring is expected to be the most perturbed 
of all the Uranian rings in this framework. The next 
largest expected amplitudes are the Cordelia 24:25 OLR 
and the Ophelia 14:13 ILR that are thought to play a 
roll in shepherding the e ring. If this is a realistic estima- 
tion of the strength of the resonances in the system, in 
the future we may be able to detect the m = —24 mode 
on the inner edge of the e ring, which was previously 
detected by French & Nicholson (1995) with occultation 
data and by Showalter (2011) with images showing the 
ring's longitudinal brightness variations. Detecting the e 
ring edge modes will first require determining the ring's 
edge positions and the removal of the larger amplitude 
m — 1 normal mode which dominates its structure. Our 
analysis of these ring residuals as well as those for the 
other rings, whose structure is dominated by previously 
known normal modes, is ongoing and will be presented 
in a subsequent publication. 


4. RESULTS 


After searching mode values from m — —25 to 25 of 
all the rings, the strongest new feature we've found is an 
m = 3 structure of the 7 ring consistent with the expec- 
tations discussed above. In Figure 4 we show the shallow 
minimum in RMS for our 7 ring m = 3 fits. The top plot 
shows the RMS deviations of the model radii from the 
observed radii at each pattern speed for m = 3, zoomed 
in on the minimum. Listed are the best fitting pattern 
speed, the semi-major axis of the Cressida 3:2 ILR, and 
the expected pattern speed for an m — 3 normal mode 
marked by the dashed line. Note that the best fitting 
pattern speed and the expected pattern speed for the 
semi-major axis of the 7 ring are offset because this is 
not a normal mode oscillation, but is instead the effect 
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Figure 3. The left mosaic shows the relative forced amplitude (darker — larger amplitude) of all first-order resonances of inner 
moons out to Perdita interacting with the Uranian rings assuming the moons each have a density of 1.3 g/cm. The largest 
amplitude resonances are labeled, while the fainter patches in the mosaic are evidence of resonances within the system that do 
not fall close enough to any rings and would have much smaller amplitudes. The actual values of the significant resonances are 
plotted on the right side to compare with the mosaic, but with the addition of a range of moon densities (0.5 to 1.3 g/cm?), 


calculated using Equation 5. 


of a resonance with a satellite whose perturbations force 
the pattern speed to match the satellite's mean motion. 
We further refine our best fit solution and formal er- 
rors by applying the best fit parameters (a, Am, Om, 
and Qp) as a set of starting parameters for MPFIT, a 
non-linear least squares fitting IDL function (Markwardt 
2009). We've initially assumed an uncertainty of 1 km 
in each of the 49 observed radii of the 7 ring, but found 
a reduced chi-squared of 0.308 << 1. We fit again to 
obtain the listed errors using the standard deviation per 
degree of freedom (c /v) as a rescaled uncertainty in our 
observed radii which better represents the error of these 
data. The bottom plot shows the best fitting model ra- 
dius curve on top of the observed radial separations from 
the fit semi-major axis of the ring, Ar = r — a. We've 
listed the final fit parameters and chi-squared analysis 
in Table 1. 

'The best fitting pattern speed for this mode, 776.58208 
+0.00169 ° /day, is strikingly close to the published mean 
motion of Cressida, the fourth moon from Uranus. Most 
recently Showalter & Lissauer (2006) listed Cressida’s 
mean motion as 776.582789 + 0.000059° /day. All three 
of the measurements of Cressida’s mean motion listed in 
Table 1 are well within the uncertainty of our detected 
pattern speed, supporting the proposed connection be- 
tween this m = 3 structure of the 7 ring and Cressida. 


To solidify that the m = 3 structure is real and is a re- 
sult of perturbations from Cressida, we have inspected 
the alignment of the structure with Cressida. In this 
case, the 7 ring (a = 47176.447) is located exterior to 
the resonance (ares = 47171.51), and the dynamical 
model predicts that one of the three outer radial ex- 
tents should track the motion of Cressida. That is, as 
the m = 3 structure and Cressida both rotate around 
Uranus at NCressida ^ Lin one of the apoapses is con- 
stantly aligned with Cressida. This can be confirmed 
by noting that the m — 3 structure has a phase offset 
63 = 58.81 + 6.12? (this is the longitude of one of the 
3 periapsis), which is roughly 60? offset from Cressida's 
longitude (359.50?) at the epoch of the fit. We show this 
alignment more precisely in Figure 5, where we have de- 
termined the offset of each occultation scan longitude 
relative to Cressida's longitude at the observation time, 
|m|(A — Acressida). The apoapse of the phase-wrapped 
structure lags the longitude of Cressida by only 6 + 11? 
(Cressida's longitude is 0? and the fit sinusoid's largest 
radial excursion occurs at 354°). This suggests that the 
perturbations on the 7 ring are due to its proximity to 
the 3:2 ILR with Cressida. 

Perhaps the most significant result of this work, shown 
in Table 2, is a determination of Cressida's mass us- 
ing Equation 5. Given As = 0.667 + 0.113 km we find 
MCressida = 2.52: 0.4 x 1017 kg. We use the effective ra- 


Table 1. 7 ring m = 3 best fit 


Parameter Final fit and scaled errors 
à (km) 47176.447 + 0.086 
As (km) 0.667 + 0.113 
ós (°) 58.81 + 6.12 
Qp (° /day) 776.58208 + 0.00169 
NCressida (° / day) 776.582789 + 0.000059" 
NCressida (°/day) 776.582414 + 0.000022” 
NCressida (^/day) 776.582447 + 0.000022° 
x? 13.861 
x? u 0.308 
a/v (km) 0.555 
N 49 
# of parameters 4 


NoTE—Listed on top are the four fit parameters 
and their formal 1-ø errors resulting from our fi- 
nal fit, where we have assumed an error of 0.555 
km for each of the observed radii of the 7 ring. 
We also list three published mean motions of 
Cressida for comparison with our pattern speed. 
'The chi-squared and reduced chi-squared below 
are from the initial fit assuming an error of 1 km 
for each radii. The unscaled errors of the pa- 
rameters in the initial fit are roughly double the 
scaled errors from the final fit, in which we have 
used the standard deviation per degree of free- 
dom as a universal error in the observed radii. 
The degrees of freedom v = N — # of fit param- 
eters. 

“From Showalter & Lissauer (2006) 


b From Jacobson (1998) 
^ From Pascu et al. (1998) 


dius for Cressida of 41 + 2 km from Karkoschka (2001a) 
to calculate a density of 0.86 + 0.16 g/cm? for Cressida. 

For our purposes, the 7 ring is outside the width of 
Cressida’s 3:2 ILR and the resulting estimation of A,, is 
reasonable, but we note that this is not necessarily the 
case for all of the other rings and resonances. Curious 
readers should note, to test the applicability of Equa- 
tion 5, we’ve calculated a resonance half-width of ~ 3.5 
km for Cressida’s 3:2 ILR using Equation 10.23 from 
Murray & Dermott (1999) along with our newly deter- 
mined mass of Cressida. The other relevant variable 
inputs can be found in Tables 1 and 2. This half-width 
is less than the 5 km separation of the resonance and 
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Figure 4. The upper plot shows the RMS deviations of the 
observed radii, r, and the model fit (Equation 2) for a range 
of pattern speeds, Qp. Listed in the top plot are the best 
fitting pattern speed and the corresponding radial location 
of the resonance, ares. The blue line refers to the expected 
pattern speed for an m — 3 normal mode oscillation. The 
lower plot shows the best fitting model (red line) and ob- 
served radii plotted vs. m0 = m[A — OQ, (t — to) — om] after 
subtracting the semi-major axis of the ring (Ar — r — a). 


0 90 180 270 360 
[IM|(A-Acressida)] mod 360 (deg) 

Figure 5. This plot shows the 7 ring structure in a ref- 
erence frame tied to Cressida. One of the three outermost 
radial extents actually tracks Cressida, the others are located 
~ 120° apart. We obtained longitudes of Cressida at vari- 
ous times using the ura091.bsp and urall2.bsp SPICE ker- 
nels, available at https://naif.jpl.nasa.gov/pub/naif/ 
generic kernels/spk/satellites/. 


ring, confirming we are justified in using Equation 5. 


Table 2. Mass and Density of Cressida 


As (km) Radius (km) a (km) 


Gres (km) 


MCressida (kg) PCressida (g cm 3) 


0.667 + 0.113 AER 


47176.447 + 0.086 47171.51 +0.03 2.5+0.4 x 1017 0.86 + 0.16 


NoTE—We list the variables needed to solve for the mass of Cressida using Equation 5. For the calculation of moressida 
we used GMuranus = 5793951.3 + 4.4 km? s7? from Jacobson (2014) and G = 6.67408 + 31 x 107H qm? kg | s^? from 
http://physics.nist.gov/cgi-bin/cuu/Value?bg. Also note 20 Mal c 1.545 when j — m — 3 for the case of the Cressida 3:2 
ILR. The listed radius needed to calculate the density of Cressida comes from Voyager 2 photometry (Karkoschka 20012). 


Note that the — 1 km resonance half-width quoted in 
the introduction was estimated by Porco & Goldreich 
(1987) and results from an approximation of the reso- 
nance half-width equation as well as a different satellite 
mass. 


5. DISCUSSION 


Since the Voyager 2 flyby of Uranus in 1986, several 
dynamicists have explored the stability of the inner Ura- 
nian moons. The moons Bianca, Cressida, Desdemona, 
Juliet, Portia, Rosalind, Cupid, Belinda, and Perdita are 
members of the most tightly packed system of moons in 
our solar system. Nicknamed the ‘Portia group’ for their 
largest member, these satellites are thought to be unsta- 
ble on short timescales compared to the age of the solar 
system. The stability of the Portia group is known to be 
highly sensitive to the masses of the individual satellites 
(French et al. 2015), which are not well constrained. In 
fact, the mass we provide for Cressida is the first direct 
measurement of an inner Uranian satellite’s mass. Past 
simulations (Duncan & Lissauer 1997; French & Showal- 
ter 2012; Quillen & French 2014; French et al. 2015) 
have relied on treating a range of possible masses for the 
inner Uranian satellites and suggest that Cressida will 
cross orbits with Desdemona in under 10° years (French 
& Showalter 2012), given our mass density. Incorpora- 
tion of our mass for Cressida should further constrain 
the timescale of satellite orbit crossing (collisions) and 
allow a future work to determine the masses of some 
of the other satellites through their resonant interac- 
tions. Strictly speaking, our density measurement does 
not necessarily represent a common density of the in- 
ner moons. However, a lower average satellite density 
will generally result in collisions occurring in the more 
distant future. 

Karkoschka (2001b) and Dumas et al. (2003) detected 
a possible water ice absorption feature in Hubble Space 
Telescope near-infrared photometry of the largest in- 
ner moon Puck. Combining this with the previously 
mentioned size estimates has formed the presumption 
that Cressida and the other inner Uranian moons are 
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Figure 6. The known densities of selected satellites in the 
solar system are plotted versus their radii. We represent the 
individual moons associated with particular planets using 
the point styles and colors labeled in the lower right of the 
plot. 


likely composed of mostly water ice with at least a ve- 
neer or contamination of dark material to explain their 
low albedo and flat gray spectra. The range in densi- 
ties of the larger Uranian moons, determined from mass 
(Jacobson et al. 1992) and radius (Thomas 1988) mea- 
surements, have provided a presumed upper limit on 
the densities of the inner moons, usually with reference 
to the least dense major moon Miranda (1.214 + 0.109 
g/cm?)!. In Figure 6 we plot our average density of 
Cressida versus radius along with other satellites in our 
solar system, after Hussmann et al. (2006). Cressida 
is about 5096 denser than the inner icy moons of Sat- 
urn with comparable radii. It may be that Cressida, 
and the Uranian rings/moons in general, have either a 


l https://ssd.jpl.nasa.gov/?sat, phys. par 


lower porosity than these Saturnian analog or they have 
higher amounts of non-icy contaminants, as inferred by 
Tiscareno et al. (2013). The contamination of denser 
and darker material may not be as high as previously 
expected, but it is substantial regardless. 

This analysis shows that there is stil information 
about Uranus’ rings and moons found in historical and 
ground based data. Still, the best means of obtaining 
the Uranian moon masses and compositions, determin- 
ing the ultimate fate of the Portia group, and under- 
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standing the intricate structure of the rings is of course 
a Uranus orbiter mission. 


We would like to thank Phil Nicholson for his in- 
sights regarding ring occultation observations and both 
Phil Nicholson and Pierre-Yves Longaretti for several 
fruitful discussions concerning the forced radial am- 
plitudes of ring particles orbiting near Lindblad reso- 
nances. We'd also like to thank our anonymous reviewer 
for helpful suggestions and comments, ultimately im- 
proving the clarity of this work. This work was sup- 
ported by the NASA Solar System Workings program 
grant NNX15AH45G. 


APPENDIX 


Included below are the occultation observation data we used in this analysis of the 7 ring. The precise numbers 
for the ring's position are derived from an analysis of the entire Uranian ring data set, including re-determined 
pole position (Pole right ascention — 77.3105814? and declination — 15.1697826?), standard gravitational parameter 
(GM = 5.793956433 x 10° km?s^?), gravitational harmonics (J> = 3.340656 x 107? and J4 = —3.148536 x 1075), and 
time offsets. The numbers therefore can deviate slightly from previously published values. We list each observation 
ID, observing location, ring plane intercept time of the relevant electromagnetic wave observed, detected mid-radius 
of the 7 ring, m = 3 fit residuals, longitude of the observation, longitude of Cressida at this time, and reference to 
publications including the observation. Longitudes are measured in the prograde direction from the ascending node of 
Uranus’ equator on the Earth's equator of the J2000 epoch. French et al. (1988) have included all observations from 
1977-1986 in their most recent fit, but more recent observations are unpublished. 
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